Ultrasonic cavity solitons 
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We report on a new type of localized structure, an ultrasonic cavity soliton, supported by large 
aspect-ratio acoustic resonators containing viscous media. The spatio-temporal dynamics of this 
system is analyzed on the basis of a generalized Swift-Hohenberg equation, derived from the mi- 
croscopic equations under conditions close to nascent bistability. These states of the acoustic and 
thermal fields are robust structures, existing whenever a spatially uniform solution and a periodic 
pattern coexist. An analytical solution for the ultrasonic cavity soliton is also presented. 
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Many systems in nature, when driven far from equi- 
librium, can self-organize giving rise to a large variety 
of patterns or structures. Although studied intensively 
for most of the last century, it has only been during the 
past thirty years that pattern formation has emerged as 
an own branch of science One of the most relevant 
features of pattern formation is its universality: systems 
with different microscopic descriptions frequently exhibit 
similar patterns on a macroscopic level. This universal 
character of pattern formation is evidenced when the mi- 
croscopic models are reduced, under given assumptions, 
to simpler equations describing the evolution of a single 
variable, the so called order parameter [2J. Order pa- 
rameter equations are usually obtained near some critical 
point, and belong to few and well known classes, such as 
the Ginzburg-Landau or the Swift-Hohenberg equations 
and their variants. They are often based on system sym- 
metries and are independent of the microscopic differ- 
ences among systems, providing a theoretical framework 
to understand the origins of non-equilibrium pattern for- 
mation 

This approach to pattern formation has been exten- 
sively applied to nonlinear optical cavities [H, 0, Hj] , such 
as lasers, optical parametric oscillators or Kerr (cubically 
nonlinear) resonators, where light in the transverse plane 
of the cavity has been shown to develop patterns with 
different symmetries (rolls, hexagons, and also quasi- 
patterns) as well as cavity solitons (CS). The latter corre- 
spond to localized solutions often resulting from a bista- 
bility between two stationary, spatially extended states 
of the system, in the form of either self-trapped switching 
waves between two homogeneous states, or isolated pat- 
tern elements embedded in a homogeneous background. 
Optical CSs are particularly interesting objects because 
of their potential use as memory bits in optical infor- 
mation processing systems [6(. However, despite the ex- 
isting analogies between optics and acoustics, acoustical 
resonators in the nonlinear regime have been much less 
explored than their optical counterparts. Furthermore, 
pattern formation studies in acoustics are almost lacking. 
The main reason lies in the weak dispersion of sound in 



common homogeneous media, which is responsible for the 
growth of higher harmonics, leading to wave distortion 
and shock formation. These effects are absent in optics, 
which is dispersive in nature. However, in some special 
cases it is still possible to avoid the nonlinear distortion 
and recover the analogies [7[. It is, for example, the case 
of sound beams propagating in viscous media character- 
ized by a strong absorption (e.g. glycerine), where sound 
velocity depends of fluid temperature, resulting in an ad- 
ditional nonlinearity mechanism of thermal origin. For 
the majority of fluids, temperature variations induced by 
an intense acoustic field result in a decrease of sound ve- 
locity, leading to a self-focusing of the beam. In the case 
of viscous fluids the characteristic length of self-focusing 
effects is much shorter than the corresponding to the 
developement of shock waves Also, high frequency 
components are strongly absorbed in such media, so in 
practice the use of a quasimonochromatic (optical) de- 
scription for wave propagation is justified. 

In a viscous medium sound propagates with a speed 
c that depends signifficantly on temperature, c = 
c (1 — <jT'), where cq is the speed of sound at some equi- 
librium (ambient) temperature, T' denotes the variation 
of the medium temperature from that equilibrium due 
to the intense acoustic wave, and a is the parameter 
of thermal nonlinearity. The propagation of sound in 
such a medium has been shown [8|, [9| to be described 
in terms of two coupled equations for pressure, p', and 
temperature, T", deviations. These equations have been 
used to address problems such as self-focusing and self- 
transparency of sound 0] . They have been also the ba- 
sis for the analysis of temporal dynamic phenomena in 
acoustic resonators [loj- In this case, a viscous fluid 
is bounded by two flat and parallel reflecting surfaces. 
One of the surfaces, vibrating at a frequency /, is an 
ultrasonic source providing the external forcing. Previ- 
ous studies on this system [lj| EH have reported, in the 
frame of the plane-wave approximation, bistability and 
complex temporal dynamics, in good agreement with the 
corresponding experiments. In this Letter we extend the 
previous model considering the effect of sound diffraction 
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and temperature diffusion in a large aperture resonator. 
These effects, which are responsible of the spatial cou- 
pling, can play an important role when the Fresnel num- 
ber of the resonator T = I 2 /XL >> 1 (being I and L its 
transverse and longitudinal dimensions, respectively). 

The intracavity pressure field is decomposed into two 
counterpropagating traveling waves, p' = p + e J (" t-fez ' + 
p_ e i(ut+kz) _|_ c c ^ wnere i j s time, z is the axial (prop- 
agation) coordinate, to — 2nf, and k = u)/cq, whose 
complex amplitudes p± are related through their reflec- 
tions at the boundaries, and the temperature field is de- 
composed into a homogeneous component and a grating 
component, T = T h + T g e i2kz + T*e~ i2kz . All these am- 
plitudes are slowly varying functions of space and time as 
the fast (acoustical) variations are explicitly taken into 
account through the complex exponentials. Under the as- 
sumption of highly reflecting plates we can adopt a mean 
field model, where the slowly varying amplitudes do not 
depend on the axial coordinate z and p- — p+ = p, end- 
ing up with the following dimensionless equations [12j | 



Tp d T P = -P + P ia + iV 2 P + i{H + G-A)P, 
d T H = -H + DV 2 H + 2 \P\ 2 , 
d T G = -t- 1 G+ DV 2 G + \P\ 2 . 
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p, H — LutpaTh, and G 



2 P § CoCp J r> " - — -p--n, 
are new normalized variables, r = t/t n is time 
measured in units of the relaxation time t n of the temper- 
ature field homogeneous component, and r p = t p /t n and 
T g = tg/th are t ne normalized relaxation times of the in- 
tracavity pressure field and the temperature grating com- 
ponent, respectively. Their original values are given by 
tp l = cqT /2L + coao, where T -C f is the transmissivity 
of the plates and «o is the absorption coefficient, and t g = 
l/4k 2 x, where \ = k /Po c p is the coefficient of thermal 
diffusivity, p is the equilibrium density of the medium 
and k and c p are the thermal conductivity and the specific 
heat of the fluid at constant pressure, respectively. Other 
parameters are the detuning A = (uj c — lo) t p , with u c 
the cavity frequency that lies nearest to the driving fre- 
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quency u>, and Pj n = ^■ 



p in , being p in the 
injected pressure plane wave amplitude, which we take as 
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real without loss of generality. Finally V 
the transverse Laplacian operator, where the dimension- 
less transverse coordinates (x, y) are measured in units of 
the diffraction length l& = cq y/t p /2u>, and the normalized 
diffusion coefficient D = xth/l^- 

The model parameters can be estimated for a typical 
experimental situation [l(| • We consider a resonator with 
high quality plates (T = O.f), separated by L = 5cm, 
driven at a frequency / — 2 MHz, and containing glyc- 
erine at 10 °C. Under these conditions the medium pa- 



rameters are Co 
1.2 x 10 3 kgm- 3 , 



2 x 10 3 ms 1 , a 
= 4x 10 3 J kg 



10 m - 



l K-\ 



Po 



and k — 0.5Wm _1 K _1 (x = 10 _7 m 2 s~ 1 ). In this case 
t p = 2 x 10~ 5 s, t e = 6 x 10~ 2 s, and our length unit 
?d = 2 mm. For a resonator with a large Fresnel num- 
ber, the relaxation of the homogeneous component of 
the temperature is mainly due to the heat flux through 
the boundaries, and can be estimated from the Newton's 
cooling law as t n ~ 10 1 s. (Remind that this is our time 
unit.) Then the diffusion constant D ~ 10°, and the nor- 
malized decay times r p ~ 10~ 6 , and r g ~ 10 -2 under 
usual conditions. We see that the problem is typically 
very stiff: < r p <C r g <C 1. In the following the re- 
sults will be given to the lowest nontrivial order in these 
smallest decay times in order to not overburdening the 
expressions. 

The spatially uniform steady state can be obtained 
by neglecting the derivatives in Eqs. |l|. Introduc- 
ing the notation W = \P\ 2 , W- m = |-Pm| 2 , one has 
G = t s W, H = 2W, and 



Wi r 



W + (A - 2Wf W. 



(2) 



The characteristic curve W vs. W nl can display an 
S-shape, typical of the optical bistability of coherently 
driven optical Kerr cavities, as we show next. The two 
turning points of the characteristic (W- m ,±i W±), verify- 
ing dWin/dVF = 0, and its inflection point (Wi n ,i, W\), 
defined by d 2 W / i n /dH /2 = 0, are given by W± — 



2A±\/A 2 -3 



and W\ 



A 
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( . , , and the corresponding val- 

ues of the input pressure follow from Eq. (J5J). Note 
that bistability requires A > Ao = \/3, the nascent 
bistability (NB) occurring at A = Ao, in which case 
W± = W\ = Wo = -7=. It is very easy to show that 



and that 



this NB inflection point occurs at Pi n ,o = 337 
the values of the model variables become Po = > 
Po = 2 Wo, and Go = t s Wq. Figure 1 shows the be- 
haviour of the homogeneous steady state illustratin g th e 
effect of detuning on the character of the solutions [13( . 

Inspired by previous works on passive optical cavities 
3 EH I we shall concentrate our study on conditions close 
to nascent bistability (A ~ Ao) and to the vicinity of 
the inflection point, (P, H, G, P n ) — (Po, Ho, Go, Pn,o), 
where interesting pattern formation properties can be en- 
visaged. Let us just comment before considering in detail 
the spatiotemporal dynamics of the thermoacoustic res- 
onator that, by focusing on the above conditions we can 
describe in a consistent way the S-shape of the charac- 
teristic curve by using standard multiple scale methods 

We consider next the stability of the spatially uni- 
form steady state. We consider perturbations of the form 
exp (ik • r + Xt) and linearize the model Eqs. {l} with re- 
spect to them. Further considering the relevant case < 
r p < r g < 1, see typical values above, one obtains that, 
to the leading order, the Lyapunov exponent A can take 
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on the following values: A G {— r p 
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FIG. 1: Characteristic curve output intensity vs input inten- 
sity as following from Eq. ((2)1 , for three values of the detuning 
A. Monostable regime (dotted line) for A = 1, nascent Insta- 
bility (full line) for A = v^, and bistable regime (dashed line) 
for A = 2 are shown. The symbol marks the position of the 
inflection point (Xo,Yo) at NB. 



where 



A+ = -1 - Dk 



4TT- 



A + fc 2 



2W 



l + (A + k 2 -2Wy 



(3) 



k = |k|, and we remind that W = \P\ , which fixes the in- 
jection pressure through Eq. ([2|). The obtention of an or- 
der parameter equation through a multiple scales analysis 
requires that a small eigenvalue exists at small wavenum- 



ber k. In our case Aj 



reads A + = a 

b- 



bk 2 



is that eigenvalue 
a = AW 



O (fc 4 ) 
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. For small k it 

A-2W _ i 
1 + (A-2W) 2 ' 



. .. 2 — D. As we will be con- 

[l+(A-2ff) 2 ] 1+(A-2W) J 

cerned with the neighbourhood of the inflection point 
under NB conditions, let us particularize our analysis to 
that point. In such case W = Wo, A = Ao and we obtain 
a — and b — D a — D, D a — Hence, A + is actually 
small for small k in the region of interest. Finally, as 
we wish to capture pattern formation asymptotically, we 
must impose that A+ can change its sign for small k. This 
imposes that the diffusion coefficient d = D — Dq ~ 0(e), 
being e a small parameter. Additionally, a careful anal- 
ysis of the homogeneous steady state reveals that, in or- 
der to describe consistently the vicinity of the NB point, 
one must consider the scalings 6 = A — A ~ 0(e 2 ) 
and n = Pi »~ J /?' 1 ~ C(e 3 ), with P inJ = ^ (l + |). 
Finally the multiple scale analysis can be performed as- 
suming p = P - P ~ 0(e), h = H - H - 0{e) and 
g = G — Gq ~ O(e), together with the slow space and 
time scales r ~ Q(e 
dard procedures [la 
as a solvability condition at the third order in e, which 
reads 



2 ) and V ~ O(e). Following stan- 
a closed equation for h is obtained 



d T h = /i 
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lhV 2 h 
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dV 2 h 



(4) 
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h to the leading order. Equation ([4]) is 
a main result of this Letter as it describes the spatiotem- 
poral dynamics of the thermoacoustic resonator in a very 
simple way. Moreover this equation can be considered as 
" universal" for our system under typical operating condi- 
tions as the numerical values of the coefficients appearing 
in it are, in the considered limit r g — > 0, very weakly de- 
pendent on r g (details will be given elsewhere). Equation 
Q is a modified Swift-Hohenberg equation, and has been 
previously introduced in the optical context for semicon- 
ductor laser s 1171 and liquid crystal light valve (LCLV) 
systems 3, 19 1. In this context the analysis of Eq. ((4]) 
has demonstrated different pattern formation scenarios, 
including front propagation, pattern competition and lo- 
calized structures. Experimental evidence has been also 
reported in [lj^ . Despite the fundamental differences be- 
tween the thermoacoustic resonator considered here and 
the referred optical systems, the derivation of a common 
order parameter equation allows to extend the optical 
predictions to the acoustic case. 

The homogeneous steady solution ho of Eq. ([J) is im- 
plicitly given by fx = — ^Sho, being multivalued 
when S > (A > Ao), in agreement with the analysis 
of the microscopic model. In the following we focus on 
the bistable region, extending between the turning points 
ho,± = ±30r- A linear stability analysis shows that, 
for d > 0, only the upper branch can be modulationally 
unstable 1181 1. giving rise to a pattern with wavenumbcr 
These results are confirmed by the numeri- 
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cal integration of Eq. ([4]), performed in a one-dimensional 
case (V = d 2 ) which describes, e.g. a physical situation 
in which the resonator has a slab waveguide configura- 
tion confining the sound in the y direction. Figure 2(a) 
shows the spatiotemporal evolution (from bottom to top) 
of an initially homogeneous distribution, where the de- 
velopement of a modulational instability is observed, in 
agreement with the previous analysis. 

In most of the multivalued domain, the extended pat- 
terns emerging at the modulational instability are un- 
stable, being a transient state. As shown in Fig. 2(a), 
neighbour maxima collide and merge, the long-term evo- 
lution resulting in a number of localized structures or 
CSs. In Fig. 2(b) the background and peak values of 
the CSs (see a profile in the inset) depending on the con- 
trol parameter \x are shown with dots. The background 
amplitude corresponds to the stable uniform state, while 
the peak value to that of the underlying extended pat- 
tern which, as stated in the introduction, is a signature 
of CSs. Note that CSs are stable and robust structures 
whenever a pattern and a uniform state coexist. 

Figure 2 reveals some important facts. Firstly, local- 
ized structures can form spontaneously from the homoge- 
neous steady state, differently from other systems which 
require a hard local excitation of the medium. Such cir- 
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FIG. 2: (a) Numerical simulation of Eq. © for S = 3 10" 3 . 
(a) The modulational instability and the further collapse re- 
sulting in an array of CSs is shown for n — 3 1CP 5 . Time 
runs from bottom to top. (b) Maximum amplitude and back- 
ground values of the CSs (dots), together with the homoge- 
neous steady solution (full line) depending on the control pa- 
rameter 11. The inset shows the profile of a CS for parameters 
in Fig. 2(a). All variables are dimensionless. 

cumstance is very relevant in the particular case con- 
sidered here, where the nature of the driving source (a 
plane, vibrating rigid surface) makes difficult to imple- 
ment experimentally a localized excitation. Certainly, 
the proposal of techniques to address CS in the acoustic 
resonator will be most relevant. 

We have also obtained an analytical, steady CS of Eq. 
@, 

h(x) — ho + h s sech 2 (x / x s ) , (5) 

where the background ho corresponds to the stable ho- 
mogeneous solution of lower amplitude. The expressions 
of the peak amplitude h s and width x s in ([5]) are quite 
involved but can be readily obtained substituting Eq. ([5]) 
into Eq. (0J . The solution given by Eq. §5§ corresponds to 
an ultrasonic CS but also applies to optical CSs in those 
systems described by Eq. (j4|). It is a singular solution, 
existing for a special value of the injected amplitude \x de- 
termined by the rest of parameters (d, 8), which is seen 
to be a particular case of the (non-analytic) CSs family 
found in the numerical study and illustrated in Fig. 2(b). 

Concluding, we have studied the spatiotemporal dy- 
namics of ultrasound in a resonator containing a viscous 
medium -a thermoacoustic resonator. The microscopic 
equations have been reduced, close to the nascent bista- 
bility point, to a single order parameter equation pre- 
viously obtained in other contexts. The analysis of the 
reduced model demonstrates the existence of bistability 



and modulational instabilities. As a consequence, the 
system is shown to support cavity solitons, correspond- 
ing to states where ultrasound is highly localized in the 
transverse plane of the resonator. Analytical solutions 
of such ultrasonic CSs have been obtained. The results 
reported here will also find application to those systems 
described by the generalized Swift-Hohcnberg equation 
which is a generic model for pattern formation. 
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